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Abstract 

The problem of critical velocities in superfluids, that is the comprehension of super- 
fluidity breakdown by flow, has been long standing. One difficulty stems from the 
existence of several breakdown mechanisms. A major advance has come from the 
observation of single 27r phase slips, which arise from the nucleation of quantised 
vortices, that is, their creation ex nihilo. The statistical properties of the nucleation 
process in both the thermal regime and the quantum regime are identified and anal- 
ysed: vortex nucleation provides a well-documented case of macroscopic quantum 
tunnelling (MQT). In particular, a close scrutiny of the experimental data obtained 
on ultra-pure ^He reveals the influence of damping on tunnelling, a rare occurrence 
where the effect of the environment on MQT can be studied. 

La Vitesse critique dans les superfluides, c'est a dire la destruction de la superflu- 
idite par I'ecoulement du fluide, pose un probleme qui perdure. Une des difficultes 
reside dans I'existence de plusieurs mecanismes pour cette destruction. L'observation 
de sauts de phase individuels de 27r, qui proviennent de la nucleation de tourbillons 
quantifies, a constitue une avancee importante. L 'identification et I'analyse des pro- 
prietes stochastiques du processus de nucleation, tant dans le regime classique que 
quantique, ont conduit a I'etude tres circonstanciee d'un cas specifique d'effet tunnel 
macroscopique. En particulier, I'examen fouille des donnees experiment ales obtenues 
avec I'helium ultra-pur a revele I'influence de la dissipation sur I'effet tunnel, donnant 
par la un exemple rare d'interaction d'un processus tunnel avec son environnement 
macroscopique. 



The critical velocity in a superfluid is the threshold above which the flow of 
the superfluid component becomes dissipative, that is, the property of super- 
fluidity is lost. This rather broad definition encompasses a number of different 
physical situations. The following overview starts with a brief description of 
the different brands of velocities that comply with this definition. It then joins 
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Fig. 1. Dispersion curve of the elementary excitations in superfluid ^He 

the main trend of this Dossier by focusing on that which involves a nucleation 
phenomenon, namely, the nucleation of superfluid vortices. 

Neither the problem of critical velocities in superfluids nor that of the nucle- 
ation of vortices are new. The former is as old as the discovery of superfluidity 
(see the monograph by Wilks [ij). The latter, flrst discussed by Vinen in the 
early sixties [2], has met a more tortuous fate. It was first thought, still is 
in some quarters, to be impossible [2j on the grounds that such an extended 
hydrodynamical object as a vortex with a finite (quantised) circulation, involv- 
ing the motion of a large number of helium atoms, would have a vanishingly 
small probability of occurring spontaneously. For classical ideal fiuids, this is 
the essence of the Kelvin-Helmholtz theorem, which states that vorticity is 
conserved for isentropic motion of inviscid fluids. More recent experiments, 
probing superflow on a finer scale of length [3ll4l|5] . have shown otherwise. 

I give below a short account of these problems, which is hardly more than 
a guided tour of the four references P[7f8|9] written by the author and his 
colleagues over the course of many years. More extended discussions can be 
found in these articles as well as more complete bibliographies. Also, a more 
comprehensive review is in preparation. 



1 Critical velocities in superfluids 

1.1 The Landau criterion 

Landau [T] explained the superfluidity of helium-4 by the sharpness of the 
dispersion curve for elementary excitations, phonons and rotons, shown in 
Fig.[T], which is a property associated with the existence of a Bose-Einstein 
condensate. Elementary excitation energy levels e(p) being well-defined, that 
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is, having a negligible spread in energy, very low-lying states, energy-wise and 
momentum-wise, are extremely scarce. An impurity, or a solid obstacle, can 
only exchange an energy e{p) at momentum p that exactly matches the energy 
of an elementary excitation of the fluid. Unless this condition can be precisely 
met, there is no dissipative interaction between the fluid and its surroundings: 
the flow is viscousless at small flow velocities. 

If the superfluid moves at velocity Vs, the energy of elementary excitations 
in the frame of reference at rest becomes e + Vg ■ p [TfTO] . The same holds 
for a moving obstacle, by Galilean invariance. If this energy turns negative, 
elementary excitations proliferate and superfluidity is lost. The condition on 
the superfluid velocity for this to happen reads: 
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The minimum value of e/p for helium is very close to the roton minimum, 
as shown in Fig.[Tl In '^He at low pressure, — 60 m/s. The Landau critical 
velocity is smaller than the sound velocity c = 220 m/s but larger than most 
critical velocities measured in various experiments. For the much less dense 
Bose-Einstein Condensed gases, which do not exhibit roton-like features, the 
minimum is the sound velocity, c = e(p)/p|p=o- 



1.2 Feynman's approach 



Feynman [T], following Onsager, realised that, not only would vorticity be 
quantised in ^He in units of the quantum of circulation = 2'Kh/mi ~ 10~^ 
cm^/s, 7714 being the mass of the helium-4 atom (which is also a property 
associated with the existence of a Bose-Einstein condensate) but that these 
vortices would be responsible for the onset of dissipation and for a critical 
velocity in the superfluid. The basic reason for this, as spelled out clearly by 
Anderson [11], is that vortices can exchange energy with the potential super- 
flow and carry this energy away, thus causing an energy loss to the superflow. 

In order to evaluate a characteristic velocity associated with this process, let us 
consider a vortex ring of radius R. Its energy E-^ and impulse Pr are expressed 
by [TMMi] 

Pr = -Kp^KiR^ , (3) 

where Oq is the vortex core radius, taken here as the superfluid coherence 
length. 
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Fig. 2. Schematic drawing of the flexible-diaphragm Helmholtz resonator. 

Let us treat such a vortex as an elementary excitation of the superfluid, which 
it rightfully is, and apply Landau's criterion. The limiting velocity is reached 
for a radius R such that E-p^/ is at a minimum, which occurs when R is 
as large as feasible, that is of the order of the channel size d. This minimum 
value sets the velocity at which vortices can start to appear and defines the 
Feynman critical velocity: 



As discussed below, v-p is much closer to experimental values than the Landau 
critical velocity for rotons. Although this agreement is heartening, it also raises 
fresh questions: how do these vortices come about? 

1.3 The phase slips 

The phase slippage experiments that were carried out starting from the mid- 
eighties PfTS] confirmed Feynman and Anderson's views on dissipation in 
superflows [11] and brought a large measure of clarification in the critical 
velocity problem [6] and in the formation of vortices in superfiuid "^He [7]. 
These experimental results and their interpretation have since been largely 
confirmed [T6|17fl8) 

Phase slips can be studied with the help of a miniature hydro-mechanical 
device, which is basically a flexible-diaphragm Helmholtz resonator as repre- 
sented schematically in Fig.O This resonator is immersed in the superfluid 
bath. The flexible diaphragm is constituted by a Kapton membrane coated 
with aluminium. In the version shown in Fig.[2l there are two openings con- 
necting the resonator chamber to the superfluid bath. One is a micro-aperture 
in which the critical velocity phenomenon takes place. The critical event con- 
sists in a sudden jump in the resonance amplitude which corresponds to an 
abrupt change in the flow velocity through the micro-aperture and a loss of 
resonator energy. These dissipation events are interpreted as resulting from 
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Fig. 3. Peak amplitude charts at 100.7 mK (top) and 201.5 mK (bottom) at saturated 
vapour pressure in ultra-pure ^He with a ^He impurity concentration less than 10~^. 
Peak absolute amplitudes during each half-cycle of the resonator motion at 12.5 Hz 
are plotted as a function of half-cycle index; time runs from left to right. The peak 
membrane amplitude is normalised to the amplitude jump of a single slip A^i. 
The traces at the very top and bottom of the graph, expanded from the main traces 
as indicated by curly arrows, show the slip sizes in signed winding numbers of the 
quantum phase (according to flow direction, in and out of the resonator chamber). 

single vortex emission, to which is associated a slip by 27i of the quantum 
phase difference across the micro-aperture, 6(p, caused by the motion of the 
vortex across the flow stream. 

The other opening is a relatively open duct and provides a parallel path to the 
superfluid along which the quantum phase remains well determined even when 
the phase slips in the micro-aperture. A quantum of circulation builds up for 
each 271 slip along the superfluid closed loop threading the two openings. The 
operation of these resonators is described in detail in the literature (see, for 
instance [5][T9ll2n]l2T]l22] ). 

The resonator is driven on resonance by an electrostatic ac-drive applied to 
the aluminium-coated flexible membrane at a constant level. In the absence of 
dissipation, the resonance motion increases linearly in amplitude under the ac- 
tion of the drive. The displacement of the membrane, which is proportional to 
the total flow in and out of the resonator, is monitored by an electrodynamic 
displacement sensor with a SQUID amplifler [23], the output of which is con- 
verted into a digital signal by a fast analogue-to-digital converter. The peak 
amplitudes of the membrane motion A^^ for each successive positive-going and 
negative-going half-cycle of the resonance are measured and stored digitally 
in real time. The raw data consist of peak amplitude charts as represented in 
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Fig.[3l Phase slips are seen as sudden drops of the peak amphtude from one 
half-cycle to the next. 

The pattern of these peak amplitude data is processed numerically to iden- 
tify all the slips and their multiplicity, and to obtain the state of quantised 
circulation trapped in the resonator loop, labelled j. The actual flow in the 
micro-aperture is the sum of the flow driven by the membrane and of the 
persistent flow threading the micro-aperture and the parallel channel, which 
depends on the quantum state of the loop j. The amplitude drop AAi caused 
by a single phase slip in a given half-cycle of the resonance corresponds to a 
change of 6(p by exactly 27i. Normalising the membrane displacement by 
Ay4i as done in Fig. [3] provides a self-calibration of the data that is indepen- 
dent of less well known quantities such as the membrane stiffness, and the 
calibration factor of the displacement sensor. 

In the following, aperture velocities are expressed by the number of turns 
S(p/27i by which the quantum-mechanical phase winds across the aperture [21]. 
The actual flow velocity averaged over the cross-section of the micro-aperture 
is proportional to 6(p, the multiplying factor being h/m^li^. The hydraulic 
length /h of the micro-aperture is of the order of 1 [im in the experiments 
shown in Fig.[3l 



1.4 Phase-slippage experimental results 

The observation of phase slips in ^He has led to a number of quite significant 
results that have shed light onto the previously indecipherable critical velocity 
problem. I summarise below the most important qualitative features and their 
implications. 

1. The critical velocity threshold, which can be seen on time charts such 
as that shown in Fig.[3l is markedly temperature-dependent down to be- 
low 200 mK and reaches a well-defined plateau below 150 mK. These 
features can been seen in Fig. [5] and will be discussed below. As the ther- 
modynamic properties of superfluid '^He are very nearly independent of 
temperature below 1 K, this observation indicates that the critical process 
in action is not governed solely by hydrodynamics. It can be suspected 
that mechanical statistics plays a leading role. 

2. Aperture-size is not found to be a relevant factor. This feature and the 
temperature dependence mentioned above are in sharp contrast with the 
Feynman critical velocity, which, according to Eq.(l4l), exhibits a well- 
characterised dependence on size and none on temperature. 

3. The actual velocity threshold for phase slips shows significant scatter 
from one slip to the next in a given sequence, as can be seen in Fig.O 
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This scatter lies much above the background noise level of detection of 
the peak amplitudes of the resonator motion. It represents a genuine 
stochastic property of the process at work, which turns out to display a 
temperature dependence similar to that of the critical velocity shown in 
Fig.l 

4. The phase slip pattern shows quite reproducible properties in the course 
of a given cool-down as long as the experimental cell is kept at a tem- 
perature below 10~15 K. If the temperature is cycled up to nitrogen 
temperature or above, small changes to the critical threshold and the 
pattern itself can occur. This is likely due to changes in the surface state 
of the cell, i.e. contamination of the micro-aperture walls by solidified 
gases. 

5. Quite importantly, phase slips are the signature that quantised vortices 
are created in aperture flow above a well-defined threshold of flow velocity. 
This statement arises from the highly reproducible phase change, which 
is measured to be very nearly 27r and to amount to changes of precisely 
one quantum of circulation in the superfluid loop threading the micro- 
aperture and the long parallel channel (see Fig.[2l). A detailed scenario for 
the occurence and development of phase slips that shows how the phase 
difference by 2tt develops has been described by Burkhart et al [25] and 
is discussed below. Different mechanims have been proposed [26ll27|28| 
for which it is unclear that the end product of the nucleation process is 
actually a vortex. 

Critical velocities and phase slips in the superfluid phases of ^He show different 
features that will be briefly touched upon in Sect IS 



1.5 Several kinds of critical velocities 

The compilation of the critical velocity data in various apertures and channels 
from various sources available in the literature presented at the Exeter Meeting 
in 1990 [fij and shown in Fig. [4] has not been updated. Two different critical 
velocity regimes appear clearly on the graph in Fig.Sl a fast regime for small 
apertures, of the phase-slip type, and a slower regime for larger channels, of the 
Feynman type. More recent data confirm this behaviour. In some occasions, 
switching between these two types of critical velocity has been observed in the 
course of the same cool-down [29f30] . 

The data points from various sources [6] for these two different types of critical 
velocity do not fall on well-defined lines as can be seen in Fig. [4] but merely 
bunch into clusters of points. As already stated, critical velocity values in 
apertures and capillaries are not very reproducible from experimemt to ex- 
periment, indicating that some less-well-controlled parameters, besides size, 
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Fig. 4. Critical velocity data vs channel width o - older data [T]; o and A - tem- 
perature-dependent and temperature-independent data, from Ref.[^. For the tem- 
perature-dependent data, the highest value, i. e. that at the lowest temperature, has 
been retained. The dash-dash line is obtained from the Feynman criterion, Eq.(l4|. 

temperature and pressure, also exert an influence. 

As a basis for comparison, it is worthwhile to also mention the findings of 
the ion propagation studies in superfluid "^He at various pressures, which have 
been reviewed by McClintock and Bowley [3TI|32| . Ions can be created in liquid 
helium and accelerated by electric fleld until they reach a critical velocity. 
The resulting drift velocities are measured by time-of-flight techniques. For 
negative ions, hollow bubbles 30 A in diameter with an electron inside, two 
different behaviours are observed: 

• Below about 10 bars, vortex rings are created, on the core of which the 
electron gets trapped: the drift velocity suddenly drops from that of the 
negatively charged bubble to the much slower vortex velocity [33] . 

• Above 10 bars, the accelerated ion runs into the roton creation barrier 
before vortex rings can be created. The Landau critical velocity is ob- 
served to decrease from about 60 m/s at SVP down to 46 m/s at 24 bars 
as the roton parameters change with pressure while the vortex creation 
velocity increases with pressure. 

• Around 10 bars, both critical velocities, the Landau critical velocity for 
the formation of rotons and that for the formation of vortex rings can be 
observed to occur simultaneously because ions can be accelerated above 
the threshold for roton emission. 

These ion propagation measurements provide a vivid illustration not only of 
the existence of a critical velocity obeying the Landau criterion but also that 
roton creation and vortex formation constitute different phenomena and can 
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exist concurrently [34]. The vortex emission threshold displays other notewor- 
thy features. It depends on temperature in a non-trivial way, comparable to 
that of the phase-slip critical velocity with the appearance of a plateau below 
~ 300mK. It also shows the marked dependence on ^He impurity concentra- 
tion observed for phase slips in micro-aperture flows but not in larger channels. 
In both ion propagation and aperture flow measurements, vortex formation 
displays very similar features. 

Altogether, a careful study of the experimental data in superfluid '^He reveals 
three different, well-deflned, types of critical velocities, one that is the cele- 
brated Landau critical velocity, another that seems related to the Feynman 
criterion with all the uncertainties on the hydrodynamical process of vortex 
creation in larger channels, and a third, for phase slips, which is in want of an 
explanation: how are the vortices of phase slips in aperture flow created, and 
how does the situation differ from that in larger channels? 

The short answer, based on qualitative evidence, is that the temperature de- 
pendence of Vc and its stochastic properties clearly point toward a process 
of nucleation by thermal activation above ~150 mK or so and by quantum 
tunnelling below. This conclusion contradicts our daily observations of the 
formation of whirlpools and eddies. It will be seen to hold in "^He because the 
nucleated vortices have nanometric size, a fact that came to be appreciated 
because of the detailed analysis of phase slippage observations that I briefly 
relate below. 



2 Phase slip critical velocity 

A more firmly established answer to the questions formulated above comes 
from a quantitative analysis of the experimental data for phase slips. These 
experiments do provide clues that, pieced together, conclusively show that, 
in small apertures, vortices are indeed nucleated by thermal activation above 
about 150 mK, and by quantum tunelling below. 

Let me begin with some preliminary remarks. A glance at Fig. [3] reveals that 
the critical velocity threshold itself needs to be defined. Also, the local value 
of the critical velocity is not measured directly. Experiments record the mean 
value of the volume flow, which is assumed to be proportional to the local 
values of the flow fleld velocity; this assumption breaks down in the presence 
of vorticity and has to be taken with a grain of salt (see Secl5]). The value 
of the critical threshold is not even reproducible from one cool-down to the 
next with the same experimental cell. This lack of reproducibility in the mea- 
surements, both in micro-apertures and in larger channels, has obscured the 
critical velocity problem for a long time. It must however be considered as an 
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integral part of the problem. 



Now, on with the real topic of this Dossier: nucleation. The first piece of 
evidence for the nucleation of vortices, that is their creation ex nihilo, rests 
on the temperature dependence of the phase-slip critical velocity shown in 
Fig. [HI which increases in a near-linear manner when the temperature decreases 
from 2 K to ~0.2 K. I mean by linear a functional dependence going as 
Vc = fo(l — T/Tq). As can be seen in Fig.[5l the data depart from this linear 
dependence below 200 mK, where they reach a plateau, and above 2 K because 
the critical velocity goes to zero at T^. 

This temperature dependence, first observed in 1985 at Orsay [5] is now a 
well-established experimental fact [8]. It came as a surprise at first because 
the critical velocities observed before were temperature-independent below ~1 
K. As the quantum fluid is nearly fully in its ground state below IK- the 
normal fluid fraction becomes less than 1 % - one is led to suspect, as was 
done in Ref.[5], that an Arrhenius-type process must come into play. If such is 
the case, that is, if a thermal fluctuation in the fluid with an energy of at most 
a few k-^T can trigger the appearance of fully-formed vortex out of nowhere, 
the energy of this vortex must also be of the order of a few kBT: it must be 
a very small vortex. But very small vortices require rather large superfluid 
velocities to sustain themselves. A careful analysis of the situation is thus in 
order. 

The nucleation rate for thermally activated process is given in terms of the 
activation energy by Arrhenius' law: 

^K^^[(l+aY''-a]exp{-^}^ (5) 

where co'o/2vr is the attempt frequency and the activation energy. The cor- 
rection for dissipation has been introduced by Kramers to describe the escape 
of a particle trapped in a potential well and interacting with a thermal bath 
in its environment. The particle undergoes Brownian motion fluctuations and 
experiences dissipation. This dissipation is characterised by a dimensionless 
coefficient a = 1/2uot, t being the time of relaxation of the system toward 
equilibrium. In superfluid helium, dissipation is small, although some dissipa- 
tion is necessary for the system to reach equilibrium with its environment. Its 
influence on the thermal activation rate is very small and will be neglected 
in the following. However, this will not be the case anymore in the quantum 
regime. 

Let us derive the expression for the critical velocity that stems from the Ar- 
rhenius rate, Eq.(l5|). In experiments performed in a Helmholtz resonator, such 
as those shown in Fig. [HI the velocity varies periodically at the resonance fre- 
quency as Vp cos(a;t), Vp being the peak velocity. The probability that a phase 
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Fig. 5. Critical velocity, normalised to the zero temperature linear extrapolation 
value vq, vs T, in kelvin: {o), Ref.[7j, for ultra-pure ^He; (a), Ref . |35| . The plain 
curves are computed from the half-ring model (see SectH]) for oq = 2.2, 3.2, 4.5, 6.0 
A and are normalised to match the experimental value at 0.5 K. The inset shows the 
influence of ^He impurities on Vc. (o), 3 ppb '^He in ^He; (a), 45 ppb, from Ref.[36j. 

slip takes place during the half-cycle uti = —njl , tutf = 7r/2 is 




- r(P, T, cos(tut') I dt' 



-2nkBT 



exp 



knT. 



(6) 



Equation ^ results from an asymptotic evaluation of the integral at the 
saddle point t = 0. The accuracy of the asymptotic evaluation ([6]) becomes 
questionable for T — > as the energy barrier vanishes. But, as we shall see, 
quantum effects take over and the energy barrier never actually vanishes. 

The critical velocity Vc is defined as the velocity for which p = 1/2. This 
definition is independent of the experimental setup, except for the occurrence 
in Eq.(l6l) of the natural frequency of the Helmholtz resonator u. The implicit 
relation between Vn and E„ then reads: 



-2nkBT 



2^0; \ dE^I dv\^ 



■ exp 



knT 



In 2 



(7) 
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Fig. 6. Statistical width of the critical velocity transition, normalised to the linear 
extrapolation limit at T = 0, uq, in terms of temperature: (o), Ref.[^; (a), Ref.[ 
(x), Ref.ETl. 



We note that, in Eq.([7]), the attempt frequency is normalised by the resonator 
drive frequency: the Brownian particle attempts to escape from the potential 
well at rate uj^j^i: but an escape event is likely only in the time window in a 
given half-cycle of the resonance during which the energy barrier stays close 
to its minimum value EjyV^^. This time interval is inversely proportional to u;, 
which explains why an instrumental parameter gets its way into Eqs.(l6|) and 



The velocity at which individual critical events take place is a stochastic quan- 
tity. Its statistical spread can be characterised by the 'width' of the probability 
distribution defined [7l|38| as the inverse of the slope of the distribution at fc, 
(9p/c}f |„J~^. This critical width is found to be expressed by: 



In 2 




dv 



+ 



1 dE, 



ksT dv 



1 -1 



(8) 



In the experiments, at low temperatures and large critical velocities, the 
quantity in curly brackets in the right hand side of Eq.((8|) is small with 
respect to the last term so that the width is simply expressed as Avc = 
-(2/ln2) ksT (^dE^/dv\^y . Thus, the statistical width is an approximate 
measure of the inverse of the slope of Eg, in terms of v. 
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Fig. 7. Left: Probability p vs slip velocity in winding number. The plain curve is a 
non-linear least square fit to the analytical form ((H), which contains two unknown 
parameters, Vc and Avc- The critical velocity is defined as that for which p = 1/2. 
Right: Nucleation rate F expressed in s~^ vs slip velocity in winding number in 
ultra-pure ^He at 17.70 mK and saturated vapour pressure on a semi-logarithmic 
scale. The line is a linear fit to the data. 



We now have precise definitions for Vc and At>c. These quantities are derived 
from p, itself obtained by integrating the histograms of the number of nu- 
cleation events ordered in velocity bins. The outcome of this procedure is 
illustrated in Fig. [3 p shows an asymmetric-S shape characteristic of the dou- 
ble exponential dependence ofp on v, Eq.^, a consequence of Arrhenius' law, 
Eq.(l5l), being plugged into a Poisson probability distribution. The observation 
of this asymmetric-S probability distribution constitutes another experimental 
clue for the existence of a nucleation process. 

The quantities Vc and Avc are easily extracted from the probability curves 
but going from fc and Afc = — (2/ln2) k^T (^dEi^/dv\^^^ back to i?a(^^) 
and u by numerical integration of the differential equation ([7]) requires more 
work and introduces additional errors. As discussed in Ref.[9j, an improved 
procedure consists in obtaining directly the escape rate F from the phase slip 
data. This quantity is the ratio, for a given velocity bin, of the number of slips 
which have occurred at that velocity to the total time spent by the system 
at that given velocity. The outcome of this procedure is illustrated in Fig. [71 
The slope of \iaF{v) directly yields dEi^/dv\^j the value of InT at Vc gives a 
combination of Ihuq and i?a(fc), which is still not easy to cleanly disentangle 

m- 

But the experiment itself offers help [40] as I now describe. 
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3 Vortex nucleation: thermal vs quantum 



Below 0.15 K, Vc ceases abruptly to vary with T, as seen in Fig. El For ultra- 
pure ^He (less than 1 part in 10^ of ^He impurities), vdT) remains flat down to 
the lowest temperatures (~ 12 mK) reached in the experiment. The crossover 
from one regime to the other is very sharp. At the same crossover temperature 
Tq, At>c also levels off sharply. It is believed on experimental grounds that this 
saturation is intrinsic and is not due to stray heating or parasitic mechanical 
vibrations; this question is of paramount importance and considerable efforts 
have been devoted to lift all uncertainties and completely elucidate the matter 

Even if all possibilities of an experimental artifact are cleared out, the mere 
observation of a plateau in Vc is no sufficient proof for a crossover from the 
thermal regime to the quantum one: the effect of ^He impurities, shown in the 
insert of Fig.[5l also gives a levelling-off of Wc vs T. This effect has been studied 
in detail in Ref.^36j and is well understood. Incidentally, it shows that the phase 
slip phenomenon taking place in the micro-aperture tracks the temperature 
down to below ~ 12 mK, the lowest temperature in these experiments: there 
is no spurious temperature saturation effect. 

If the nucleation barrier were undergoing an abrupt change at Tq, for instance 
because of a bifurcation toward a vortex instability of a different nature [26] , 
in all likelihood Avc would jump to a different value characteristic of the new 
process (presumably small since Vc reaches a plateau). Such a jump is not 
observed in Fig.O Furthermore, Vc levels off below Tq, which would imply 
through Eq.Q that Ea becomes a very steep function of v, but Avc also levels 
off, which, through Eq.^, would imply the contrary. This remark leads us to 
investigate the possibility that, below Tq, thermally-assisted escape over the 
barrier gives way to quantum tunnelling under the barrier [41]. This transition 
would induce plateaus below Tq for both Vc and At>c. 

I state again that, during the course of these investigations at Orsay-Saclay, 
the group of Peter McClintock at Lancaster concluded in their ion propagation 
studies to the existence of a crossover around 300 mK from a thermal to a 
quantum regime for the nucleation of vortices [42], as predicted by Muirhead, 
Vinen, and Donnelly [4]. There certainly are significant differences between 
the ion limiting drift velocity and aperture critical fiow - in particular, the 
latter is nearly one order of magnitude smaller - but the qualitative similarities 
are striking. We thus have two completely different types of experiments that 
point toward vortex nucleation, both in a thermal regime and in a quantum 
one. 
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3.1 The macroscopic quantum tunnelling rate 



To proceed with our investigation, let us now make the assumption that below 
Tq, zero point fluctuations do take over thermal fluctuations. The potential 
barrier is not surmounted with the assistance of a large thermal fluctuation, 
it is tunnelled under quantum-mechanically; non-conservation of energy is 
not a problem if it is brief enough, as stated by the Heisenberg uncertainty 
principle for energy. The quantum-tunnelling event is "assisted" by the zero 
point fluctuations [43]. What is remarkable here, and not necessarily easy to 
admit, is that such an energy non-conserving process does affect a macroscopic 
number of atoms, that necessary to form a vortex of about 50 A in length, as 
we shall see below. 

Such "macroscopic quantum tunnelling" (MQT) processes have been the object 
of numerous experimental and theoretical studies, mainly in superconducting 
Josephson devices. The case of vortices in helium can be treated in a very 
similar manner, as done in Ref.[9|. Before giving a brief relation of MQT for 
vortices in ^He, I summarise some of the basic results of the extended body of 
theoretical studies that followed Caldeira and Leggett's original work [44] . 

The quantum tunnelling rate of escape out of a potential well V{q) is a text- 
book problem [45]. The rate is proportional to exp —S/h, S being, in the WKB 
approximation, the action of the escaping particle along the saddle-point tra- 
jectory at the top of the potential barrier, the so-called "bounce" [46]. For 
a particle of mass m and energy E escaping from a one-dimensional barrier 
V{q), this action reads 



The determination of the bounce yields the points qi and q2 at which the 
particle enters and leaves the barrier. 

A discussion of the quantum tunnelling of vortices thus requires a Lagrangian 
formulation of vortex dynamics. Such a formulation has been carried out in 
particular by Sonin [Hj (see also Ref.[47] for an extended discussion). However, 
analytical results can be obtained only at the cost of approximations and yield 
less than fair comparison with experiments (see the discussion in Ref.[8]). 

Here, I follow, as in Ref.[9], the usual approach taken in the literature for 
Josephson devices [44f48] . which is to choose for V{q) a simple analytic form 
limited to a parabolic and cubic term in q: 




(9) 




(10) 
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where luq is the angular frequency of the lowest mode of the trapped particle 
and gb the generalised coordinate of the barrier top location. The barrier height 
-Eb is equal to mujlq^/Q. 

This simple form is of general applicability when the applied velocity is close 
to the limit, which I call fco, where the energy barrier vanishes and the system 
"runs away", the so-called "lability" point. At this point, the critical velocity 
is reached even in the absence of thermal or quantum fluctuations. Such a 
hydrodynamic instability threshold at which vortices appear spontaneously 
has been shown to occur in numerical simulations of flows past an obstacle 
using the Gross-Pitaevskii equation by Frisch et al. [49] and others [50|51f52) . 

The zero-temperature WKB tunnelling rate for the cubic-plus-parabolic po- 
tential ^^b, Eq.Jin]), is found to be [44] 



the action 5*0 being equal to 36-Eb/5ci;o- 

From this result, we may anticipate that the crossover between the quantum 
and the thermal regime lies around a temperature close to that for which the 
exponents in Eqs.(l5]) and (fTTl) are equal, namely T = Scuq/SG/cb - assuming 
that the activation energy in Eq.(l5|), -Ea, reduces to the simple cubic-plus- 
parabolic form, Ey,. A more precise study of the mathematical properties of 
the quantum channel for escape leads to the following relation [53f54] 



Once the crossover temperature has been determined from experiment, Uq 
is flxed to pinpoint accuracy compared to the fltting procedure outlined in 
the previous Section. This is where we get help from the experiment because 
both the value of uo is now completely pinned down and the interpretation 
of the experiment in terms of a nucleation process is confirmed. The values of 
the barrier height E^, at each given velocity then follow easily, using the full 
expressions for the rate in terms in terms of i?b, and, also, the damping 
parameter a. 

3.2 Friction in MQT 

Damping turns out to matter significantly in the quantum tunnelling of semi- 
macroscopic objects. The relevance and applicability of the concept of quan- 
tum tunnelling to macroscopic quantities such as the electric current through 
a Josephson junction or the fiow of superfiuid through a micro-aperture, al- 
though still sometimes questioned, have been checked in detail for the Joseph- 
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Fig. 8. Potential well trapping a particle in one dimension. The particle can escape 
to the continuum of states to the right. The lowest mode at the bottom of the well 
has angular frequency wq; Wb would be the corresponding quantity if the potential 
was inverted bottom over top. There can exist intermediate energy levels within the 
well, which are populated according to the Boltzmann factor. 



son effect case [55]. One of the conceptual difficulties, besides the large num- 
ber of particles involved, is that the macroscopic system is coupled to an 
environment that acts as a thermal bath; this coupling gives rise to a source 
of fluctuations and semi-classical friction. This issue was tackled by Caldeira 
and Leggett [44|, and a number of other authors (see, for instance, [53|54] ). 
In the case of vi^eak ohmic damping (a <C 1) and for the cubic-plus-parabolic 
potential, the tunnelling rate takes the form [44|56f57) : 



r. 



qt 
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(13) 



Thus, according to Eg. ([131) . damping depresses the MQT escape rate at T = 
- a is a positive quantity - and introduces a temperature dependence that in- 
creases the rate as T increases. These effects are large, even for weak damping, 
because they enter the exponent of the exponential factor in Eq.(fT3ll. Relation 
( fT2ll between Tq and ojq is nearly unaffected by damping: uq is simply changed 
into t^o[(l + a^Y^"^ — a] according to Eq.(l5|), a minor modiflcation for a <^ 1. 



Equation ( fTSll is valid up to about Tq/2. From Tq/2 to ~ Tq, one has to 
resort to numerical calculations [56]. In the thermal activation regime, T > Tq, 
quantum corrections affect the Kramers escape rate up to about 3Tq and 
can be evaluated analytically. These high-temperature quantum corrections 
depend only weakly on friction. A complete solution of the problem of the 
influence of friction, weak, moderate or strong, in the regime where thermal 
fluctuations still prevail but quantum corrections cannot be neglected has 
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Fig. 9. Logarithm of the escape rate normahsed to the attempt frequency in terms 
of inverse temperature, also normalised to ujq for various value of the damping pa- 
rameter, adapted from Ref.[56]. The dot-dot line is a hand-sketch of the situation 
where a increases with temperature, starting from zero at T = 0. 

first been worked out in the classical regime (T ^ Tq) by Grabert [58]. The 
extension to the temperature range T > Tq was then carried out by Rips and 
Pollak [59] who showed that the rate for arbitrary damping in the temperature 
range T > T„ can be factorised into physically meaningful terms: 



(14) 



namely, the classical Kramers rate Tk, the quantum correction factor /q, and 
the depopulation factor T . The high temperature limit of /q is 



(15) 



in which ojq and are the confining potential parameters depicted in Fig. [3 
Analytic results for /q are known to slightly below Tq [56|60| . 

The depopulation factor T describes the fact that the escape process even- 
tually depletes the occupancy of the energy levels inside the potential well. 
This factor, the expression of which is too bulky to be reproduced here (see 
Ref.[59]). is unity at large a when the coupling of the Brownian particle with 
the thermal bath is large. It decreases to zero as a — and the system be- 
comes effectively decoupled from the environment. In the quantum regime, 
dominated by zero point fiuctuations, level depletion does not take place and 
T is unity. For the nucleation of vortices, friction turns out to always be both 
sufficient and not too large so that depopulation corrections remain small and 

r ~ 1. 



The escape rate calculated for three values of the damping parameter a over 
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34 36 38 40 

V [winding number] 

Fig. 10. The barrier energy -Eb in kelvin vs v, the velocity in the aperture expressed 
in phase winding number obtained: (■), from the LT data transformed using the 
numerical tables in Ref.|56|: (o), from the high temperature data. The high T and 
low T analyses yield consistent results in the region where they overlap. 

the full temperature range is shown in Fig.[9l A hand sketch shows the influence 
of a temperature dependence in the damping coefficient. Such a situation is 
found in the nucleation of vortices in ^He as I now describe. 



3.3 Experimental energy barrier and damping coefficient 

From this knowledge of the theoretical analytical and numerical expressions for 
the rate T, obtained for the cubic-plus-parabolic potential, it becomes possible 
to extract from the measured nucleation rate and crossover temperature the 
values of the energy barrier in terms of Vc- The value of c^o given by Eq.(fT2il 
(uJo/271 = 2 X 10^° Hz for Tq=0.147 K) is consistent with the attempt frequency 
appropriate to the thermally-activated regime [39j and that found directly 
from the fits to the probability p as shown in Fig. [3 Furthermore, it agrees 
well (for ao = 4.5 A) with the eigenfrequency of the highest Kelvin mode that a 
vortex filament in "^He can sustain, uj+ = Ki/n-al = uq. The final step consists 
to extract the values for the energy barrier i?b from the measured escape rate. 
These values of i?b in the case of the experiments on ultra-pure "^He analysed 
in Ref.[9] are shown in Fig.fTOl 

The self-consistency of the procedure can be checked by using the values of ti^o 
and i?b derived from this analysis of the nucleation rate to compute Vc and Avc 
using Eqs.(l7l) and ((HI), mutatis mutandis, and compare with the experimentally 
determined values. We thus have a form of closure procedure to check the 
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1/T [K ] 



Fig. 11. lnr{v) vs 1/T, F being expressed in s~^ and T in K: (o), as measured at 
varying T and Vc] (*), corrected for the change of the velocity with T as explained 
in the text. The raw data from the run with the second sample of ultra-pure ^He, 
(a), (a), agree very well with that of the first sample. 
In the inset, In r{v, 



has been obtained with smoothed values of Vc- The curves 
represent the calculated values of In r{vq) with a = (dash-dash) or varying with T 
(plain) as explained in the text. The dot-dot curve is the extrapolation to l/T = 
of a linear fit to the high temperature portion of the data. 

analysis, from which we conclude that our assumption according to which 
vortices are nucleated by quantum tunnelling below Tq shows full consistency 
with the thermally- assisted nucleation regime that prevails above Tq. 

The quantitative analysis can be carried out one step further by construct- 
ing a Arrhenius plot from the experimental data and comparing directly the 
outcome to the results from theory. Arrhenius plots are drawn at constant i?b 
and varying temperature. Our results here are obtained at velocities that vary 
with temperature, hence at varying -Eb's. As can be noted in Fig. [HI the raw 
experimental, velocity-dependent, rates exhibit little variation over the range 
of parameters: escape rates are only observed in a certain window determined 
by experimental techniques. At low temperatures, T < Tq, the critical velocity 
is close to its zero temperature limit Vq and the corrections to T are small. As 
T increases above Tq, Vc decreases and T has to be determined by piecewise 
integration of dlnT/df. The high temperature extrapolation for T obtained 
in such a manner does display the expected 1/T dependence, as seen in Fig.fTTl 

The low temperature corrected T shows, as can be seen in the inset of Fig.fTTl 
a small, but real, drop below its zero temperature limit as the temperature is 
raised, thereby simply following the trend of the measured T. As illustrated 
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in Fig.[9l this drop reveals the influence of damping. A damping coeflicient a 
that increases from OatT = Oto~ 0.1 around Tq and more slowly above 
accounts for the observed drop [9]. This T-dependent dissipation also makes 
the crossover between the thermal and the quantum regimes even sharper than 
for a = 0, and closer to observations. The nucleation of vortices in "^He thus 
offers a rare observation of the effect of damping on MQT. 



4 The vortex half-ring model 



As described in the previous Section, there is strong evidence that the experi- 
mental features of the phase slip data result from quantised vortex nucleation. 
The nucleation barrier is of the order of a few kelvins (see Fig.[TOl). The 
attempt frequency ~ 2 x 10^° Hz is of the order of the highest Kelvin waves 
mode. In this Section, I wish to describe a simple model that will account for 
the features described above. This model, the nucleation of vortex half-rings 
at a prominent asperity on the walls, finds its roots in the work of Langer, 
Fischer, and Reppy [3ll6T] . Volovik [62j, and Muirhead, Vinen and Donnelly 
[4]. It was further developed and put on firm experimental findings in Ref.fTj. 

The model premise is quite simple. Consider, as done by Langer and Reppy 
in Ref.[3j, the homogeneous nucleation of a vortex ring in a homogeneous fiow 
Vs- When the ring has reached radius i? in a plane perpendicular to the fiow, 
its energy in the laboratory frame, where the observer is a rest and sees the 
superfiuid moving at velocity Wg, is expressed by 

E, = E^- P^v, . (16) 

The rest energy E^i and impulse Pr of the vortex ring are given by Eqs.(l2]) 
and ([3]). The minus sign in the right hand side of Eq.fflGl) arises because I have 
implicitly assumed that the vortex opposes the fiow, i.e., that its impulse Pr 
points straight against Vs'. this configuration minimises E^. 

The rest energy E^ increases with vortex size as PlnP and the impulse Pr 
as P^: the impulse term becomes dominant at large radii and causes E^ to 
become negative. The variation of E^ in terms of P has the shape of a confining 
well potential, which becomes shallower and shallower with increasing t>s, as 
depicted in Fig. [121 The barrier height can easily be computed numerically and 
substituted into the expression for fc, Eq.(l7]). An analytical approximation 
for Vc involving the neglect of logarithmic terms and valid for large vortices 
(P ^ ao) has been given by Langer and Reppy [3]. 

Such a critical velocity would be for the formation of a mist of vortices in the 
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Fig. 12. Energy barrier vs R, the vortex radius, in the simple half-ring model 
at various superfluid flow velocities. This figure is taken from Ref.HT] and applies for 
the model and with the units used in that reference. It is given here for illustration 
purposes, namely that potentials of the form of Eq. lfTHI) give rise to confining wells. 

bulk of the superfluid. However, this sort of vorticity condensation does not 
take place for two reasons. Firstly, the velocity of potential flows, which follows 
from the Laplace equation, reaches its maximum value at the boundaries, not 
in the bulk. Secondly, the nucleation of a vortex half-ring at the boundary itself 
involves a half of the energy given by Eq. (fT6l) . Hence, half-ring nucleation at 
the wall is always much more probable at the same velocity Vg than full ring 
nucleation in the bulk. Half of the energy for the half-ring holds for classical 
hydrodynamics, the other half being taken care of by the image in the plane 
boundary. For a superfluid vortex, the actual energy of a half-ring is smaller 
than in the classical ideal fluid because the superfluid density is depleted at 
the solid wall and the core radius increases. This effect strengthens the case 
for half-rings, as discussed in Ref.|25|. 

The barrier height can easily be computed and substituted into the expressions 
for Vc and Avc, Eqs.(l7l) and ([HI). Critical velocities Vc and statistical widths 
Avc computed in such a manner are shown as a function of temperature in 
Figsl5]and[6]for several values of the vortex core parameter ao- A value of 4.5 
A gives near-quantitative agreement with the experimental observations over 
the entire temperature range. This value of Oq is compatible not only with the 
temperature variations of Vc and Avc but also with the magnitude of the local 
Vc found to be 20 ~ 22 m/s using '^He impurities as a local velocity probe 
[36j. It exceeds that in the bulk (oq — 2.5 A), which is thought to reflect the 
proximity of the wall as discussed in greater details in Ref.[8]. With this value, 
the nucleating half-ring has a radius of approximately 15 A at the top of the 
barrier. 

Once nucleated, the vortex floats away, carried out by the superfluid stream at 
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Fig. 13. 2D (left) and 3D (right) views of the vortex half-ring trajectory over a 
point-like orifice in an infinite plane. The dash-dash lines on the 2D plot are the 
potential flow lines that emerge from the orifice. 

the local superfluid velocity and by its ov^n velocity, = dE^/dP^. It can be 
noted that, at the top of the barrier, dE^/dR = and the vortex self- velocity 
fR exactly balances the applied Vs'. the nucleating vortex is at a near standstill. 

If the flow is uniform, with parallel streamlines, nothing much happens; the 
vortex wanders away and the interaction with the normal fluid and with the 
wall causes a loss of vortex energy that eventually leads to its disappearance. 
If the flow is divergent, as in Fig.[T3l the vortex tends to follow the local 
streamlines and grow under the combined action of the potential flow and its 
own self-velocity: it then gains energy at the expense of the potential flow. In 
such a way, it can expand from nanometric to micrometric sizes and above. 
The vortex in its motion away from the micro-aperture takes a flnite lump of 
energy to remote places of the cell. This energy loss reduces the Helmholtz 
resonance amplitude in a fairly sudden manner. Such a dissipative event gives 
the signature of single phase slips that is seen in Fig.O 

This scenario for a phase slip involves a change of the phase difference between 
the two sides of the micro-aperture of exactly 27i because the vortex ends up 
crossing all the streamlines, as pictured in Fig. [131 This crossing causes the 
velocity circulation to change by exactly one quantum K4 on all the superfluid 
paths extending from one side of the aperture to the other. 



5 Pinning, vortex mills, collapses and all that 

Single phase slips are observed in experimental situations which may be loosely 
characterised as "clean", that is, for uncontaminated apertures of relatively 
small sizes (a few micrometres at the most), with low background of me- 
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chanical and acoustical interferences, etc and with probing techniques that 
do not manhandle the superfluid, namely, with low frequency Helmholtz res- 
onators. When these conditions are not met, flow dissipation occurs in a more 
or less erratic manner in large bursts - multiple phase slips or 'collapses' of 
the superflow. 

Multiple phase slips and collapses constitute an apparent disruption of the 
vortex nucleation mechanism described in the previous Section. Their proper- 
ties have been studied in detail in Ref.[22] and are briefly mentioned below, 
together with possible mechanisms for their formation. It is likely that these 
events provide a bridge between the "clean" single phase slip case and the 
usual situation of the Feynman type critical velocities that are temperature- 
independent below 1 K and dependent on the channel size. This problem, 
which is not fully resolved at present, almost certainly involves some form of 
preexisting vorticity. 

5.1 Remnant vorticity and vortex mills 

Remnant vorticity in ^He, which has long been assumed, has been shown 
directly to exist by Awschalom and Schwarz by looking at the trapping of 
ions by vortex lines [63]. Vortices, presumably nucleated at the X transition 
where the critical velocity is very low, remain stuck in various places of the 
superfluid sample container. This trapped vorticity, according to Adams et al. 
[64] . either is quite loosely bound to the substrate and disappears rapidly, or 
is strongly pinned and is dislodged only by strong perturbations. 

To account for laboratory observations and with the outcomes of numerous 
numerical simulations of vortex dynamics, Schwarz has proposed the following 
formula for the velocity at which vortices unpin [65] . 



D being the size of the pinned vortex and b being a characteristic size of the 
pinning asperity. Equation (fTTI) bears a strong resemblance with that for the 
Feynman critical velocity, Eq.Q. Long vortices unpin at very low velocities 
unless they are perched on a tall pedestal, but very small vortices pinned on 
microscopic defects at the cell walls can in principle exist under a wide range 
of superflow velocities; a 200 nm long vortex fllament pinned at both ends on 
20 A asperities resists transverse flows of velocities up to 10 cm/s. 

In connection with the critical velocity problem, the long standing sugges- 
tion by Glaberson and Donnelly [GGj of vortex mills still prevails. In these 
authors' views, imposing a flow on a vortex pinned between the opposite lips 




(17) 
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of an aperture would induce deformations such that the vortex would twist on 
itself, undergo self-reconnections, and mill out free vortex loops. First, we note 
that, according to Eq.lfTTl). such a mill must involve a pinned vortex of sub- 
micrometric size for any flow velocity above ~ 1 cm/s in order for the pinned 
vortex not to be washed away. Thus, it cannot as such account for Feynman- 
type critical velocities found in large channels. Also, as shown by numerical 
simulations of 3D flows involving few vortices only [67j, vortex loops and fil- 
aments are stable even against large deformations. Vortices are not prone to 
twist on themselves and foster loops. It takes the complex flow fields associ- 
ated with fully developed vortex tangles to produce small rings [68|69] . And it 
takes some quite special vortex pinning geometry to set up a mill that actually 
works. 

Schwarz has demonstrated the existence of such a mill by numerical simula- 
tions [70]. Let us take a vortex pinned at one end and fioating along the fiow 
streamlines with its other end moving freely on the wall; this vortex develops 
a helical motion, a sort of driven Kelvin wave, and reconnects sporadically to 
the wall when the amplitude of the helical motion grows large enough. This 
helical mill, which has to be of sub-micrometric size to stand the fiow, does 
churn out fresh vortices. 

Vortex mills are thus unlikely to explain the critical velocities of the Feynman 
type in the simple scheme suggested by Glaberson and Donnelly [66] • However, 
the occurrence of multiple slips, which can be seen in Fig.[3l is probably caused 
by some form of vortex mills on a microscopic size. Before coming to this topic, 
I need to describe multiple slips in greater details. But, at this point, the above 
remarks on the stability of vortex loops or half-loops in their course already 
make it unlikely that multiple slips be due to the production of small rings 
by the nucleating vortices twisting on themselves a la Glaberson-Donnelly, 
as suggested by Amar et al. [71]. We have to dig a little further to devise a 
scheme that works. 



5.2 The two types of large slips 

Besides the usual single slip pattern, there appears in Fig. [3] occasional double 
slips (i.e. involving phase changes by in) and infrequent triple slips. Raising 
the temperature to 80 mK, again for this particular cool-down, causes these 
multiple slips to occur much more frequently and to involve more circulation 
quanta on the mean. These features are described in detail in Ref . [22] . As the 
probability for a one-slip event per half-cycle is not large, that for a double slip 
is small, and it becomes negligible for higher multiples. A separate mechanism 
for their formation must be found. 
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Fig. 14. Mean size of (positive) multiple slips vs velocity in phase winding number 
in nominal purity ^He (100 ppb ^He): (a) pressure sweep from 0.4 to 24 bars at 81.5 
mK (all even values of the pressure P, and 0.4, 1, 3, 5, 7 bars) - (o) temperature 
sweep at 16 bars - (o) temperature sweep at 24 bars - (*) drive level sweep at 24 
bars, 81.5 mK - (■) temperature sweep at bar. For the temperature sweeps, from 
14 to 200 mK approximately, v first increases, reaches the quantum plateau and 
then decreases, as shown in the insert of Fig.[5l Lines connect successive data points 
in the temperature and pressure sweeps. 

Some degree of understanding of the formation of multiple slips can be gained 
by plotting the mean value of the phase slip sizes, expressed in number of 
quanta, against the flow velocity at which the slips take place [72]. This flow 
velocity is close to the critical velocity for single phase slips, i.e. the vortex 
nucleation velocity; it varies with temperature, pressure, and resonator drive 
level. A plot summarising these variations is shown in Fig. [14] for <n^>, i.e. 
in the flow direction conventionally chosen as the + direction. Slips in the 
opposite (— ) direction behave qualitatively in the same manner but the phe- 
nomenon displays a clear quantitative asymmetry. As can be seen in Fig.[T4l 
the mean slip size decreases, as does the nucleation velocity, on either side of 
the quantum plateau - a ^He impurity effect on the low T side - a thermal 
effect on the high-T side. However, it increases with pressure, contrarily to the 
nucleation velocity which decreases with increasing pressure. 

We conclude from the organisation of the data with the various parameters 
in Fig.[T4] that the magnitude of the superflow velocity does not directly con- 
trol, by itself, the occurrence of multiple slips. In turn, this implies, as will 
be discussed further below, that the phenomenon under study is not purely 
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ruled by hydrodynamics in the bulk of the fluid but involves some complex 
interplay with the boundaries. As shown in Fig. [HI the velocity threshold for 
the appearance of multiple slips depends on hydrostatic pressure; in fact, the 
P-dependence of the upturn of <n^> vs f exactly tracks that of the critical 
velocity for single phase slip nucleation. This indicates that multiple slips ap- 
pear because of an alteration, or as a consequence, of the nucleation process 
itself. 

The pattern of formation of multiple slips changes from cool-down of the cell 
from room temperature to cool-down but remains stable for each given cool- 
down. It seems to depend on the degree of contamination of the cell, degree 
which cannot easily be controlled experimentally. The detailed microscopic 
configuration of the aperture wall where nucleation takes place probably plays 
an major role in multiple slip formation. 

Another kind of very large drops in the resonance amplitude of the resonator 
was also observed, which sometimes resulted in a complete collapse of the 
resonance. These "singular" collapses, first observed by Hess j73j, occur at 
fiow velocities that are lower than the critical velocity for phase slips, some- 
times as low as t^c/S. Multiple slips are different from "singular" collapses and 
the underlying mechanisms responsible for both phenomena are bound to be 
different, as discussed below. 



5.3 In-situ contamination by atomic clusters: pinning and collapses 

In a series of experiments conducted at Saclay [741175] . in which the experi- 
mental cell was deliberately heavily contaminated by atomic clusters of air or 
H2, we observed that numerous multiple slips and collapses of the "singular" 
type occurred. The peak amplitude charts of the resonator became mostly 
impossible to interpret, except in a few instances where two apparent criti- 
cal velocities for single phase slips were observed. The higher critical velocity 
corresponds to the one observed in the absence of contamination. The lower 
critical velocity is thought to reveal the influence of a vortex pinned in the 
immediate vicinity of the nucleation site. This vortex induces a local velocity 
which adds to that of the applied flow and causes an apparent decrease in the 
critical velocity for phase slips. Because of this change, the presence of the 
pinned vortex could be monitored, the lifetime in the pinned state and the 
unpinning velocity could be measured, yielding precious information on the 
pinning process. 

This observation, reported in detail in Refs.[74] and [75], shows that pinned 
vorticity can contribute to the nucleation of new vortices at the walls of the 
experimental cell. Such pinned vortices as the one described above can, instead 
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of interacting with the nucleation site, set up a transient vortex mill of the 
helical type and generate a burst of vortices. The existence of such pinned 
vortices is established; that they can form a micro-mill is highly plausible. We 
thus have a possible explanation for multiple slip formation [8j. The pinning 
event would take place immediately after nucleation when the velocity of the 
vortex relative to the boundary is still very small and the capture by a pinning 
site easy. The micro-mill remains in activity as long as the flow is sufficient to 
maintain the helical instability, which depends on the pinning stand geometry. 
As it is set up to withstand one flow direction, it is destroyed when the flow 
velocity reverses itself in the resonance motion. It eventually re-establishes 
itself during a subsequent resonance cycle, causing a new multiple slip. This 
process depends on the precise details of the pinning site configuration and 
of the primordial vortex trajectory, factors which allow for the variableness of 
multiple slips on contamination and pressure. 

In the same experiments, we observed that a large number of unpinning events 
were taking place at an "anomalously low" unpinning velocity. A parallel can 
be made [TS] with the singular collapses that also occur at "subcritical" veloci- 
ties. In fact, both phenomena were seen quite frequently in these experiments, 
suggesting that they have a common cause. Noting furthermore that pinning 
and unpinning processes were quite frequent, releasing a fair amount of va- 
grant vorticity, it appears quite plausible that both singular collapses and low 
velocity unpinning events are caused by vagrant vortices hopping from pin- 
ning sites to pinning sites, ending up passing by close to a pinned vortex or a 
vortex nucleation site, and giving a transient boost to the local velocity, which 
pushes a pinned vortex off its perch or causes a burst of vortices to be shed. 

These observations, albeit incidental, have important consequences for the 
critical velocity problem: existing vortices, either pinned or free-moving, can 
contribute to the nucleation of new vortices at the walls of the experimental 
cell at apparent velocities much lower than the critical velocity for phase slips. 
We are thus provided with a mechanism by which superflow dissipation sets 
in at mean velocities on the large scale much smaller the velocity for vortex 
nucleation on the microscopic scale, possibly bridging the gap between phase 
slip and Feynman type critical velocities. Vortex nucleation at the walls is 
also quite likely to take part in the build-up of self-sustaining vortex tangles 
forming superfluid turbulence, up to now attributed solely to reconnection 
mechanisms 1761 . 




To conclude, the critical velocities in superfluids that are true and proven in- 
clude the Landau critical velocity for roton creation, the formation of vortices 
by a hydrodynamical instability in BEC gases [77] and in ^He [78], the nucle- 
ation of vortices by thermal activation and quantum tunnelling in '^He, both 
for ion propagation and in aperture flow. I have presented rather compelling 
experimental evidence for the interplay between vortex nucleation and pinned 
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vorticity on a microscopic scale; this evidence points toward the existence of 
helical vortex micro-mills that can generate vortices at fairly low applied ve- 
locities. Finally, vagrant vortices interacting with these mills, or with vortex 
nucleation sites, are found to generate enough vorticity to completely kill the 
applied superflow and explain singular collapses. The study of phase slippage 
has taken us a long way toward an explanation of critical velocities in super- 
fluid helium-4. 
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